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Abstract 

The existence of topological solutions and axial anomaly open a possibility 
of V and CV violation in QCD. For a reason that has not yet been estab- 
lished conclusively, this possibility is not realized in strong interactions - 
the experimental data indicate that a global V and CV violation in QCD 
is absent. Nevertheless, the fluctuations of topological charge in QCD vac- 
uum, although not observable directly, are expected to play an important 
role in the breaking of Ua(X) symmetry and in the mass spectrum and other 
properties of hadrons. Moreover, in the presence of very intense external 
electromagnetic fields topological solutions of QCD can induce local V— and 
CV— odd effects in the SU C (2>) x U em (l) gauge theory that can be observed in 
experiment directly. Here I show how these local parity-violating phenomena 
can be described by using the Maxwell-Chern-Simons, or axion, electrody- 
namics as an effective theory. Local V— and CV— violation in hot QCD 
matter can be observed in experiment through the "chiral magnetic effect" 
- the separation of electric charge along the axis of magnetic field. Very 
recently, STAR Collaboration presented an observation of the electric charge 
asymmetry with respect to reaction plane in relativistic heavy ion collisions 
at RHIC. 
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1. Introduction 

At present Quantum Chromo-Dynamics (QCD) is firmly established as 
the theory of strong interactions. At short distances the QCD running cou- 
pling constant is small due to the asymptotic freedom [1, 2]; because of this, 
the amplitudes of hard parton scattering can be reliably evaluated in per- 
turbation theory. Moreover at weak coupling one expects to find a role for 
the classical solutions of QCD equations of motion as the quantum correc- 
tions to these solutions are under control. At strong coupling, we still lack 
a reliable theoretical method for performing analytical calculations in QCD. 
Nevertheless the results of numerical lattice calculations as well as analytical 
results in related theories suggest that some of the phenomena induced by 
the presence of classical solutions persist even when the coupling constant 
grows large. 

The non-Abelian nature of QCD has very important implications for the 
structure of the classical solutions. For example, the classical Yang-Mills 
equations (non-Abelian analogs of Maxwell equations) possess non-trivial 
vacuum solutions - instantons - corresponding to the mapping of the SU (2) 
subgroup of the gauge group SU (3) onto the group of three-dimensional ro- 
tations S3 [3]. Instantons thus couple rotations in space to rotations in the 
space of color. In the presence of fermions, this property of instantons causes 
non-conservation of chirality that would otherwise be forbidden by the con- 
servation of angular momentum - but the coupling of angular momentum to 
color allows to compensate the flip of spin by a rotation in color space. In 
Euclidean space-time, instantons are static localized objects. It is well known 
that the amplitude A tun of quantum tunneling in Minkowski space is deter- 
mined in the quasi-classical approximation by the action S c i of the classical 
Euclidean solution, A tun ~ exp(—S c i/h). Likewise, in Minkowski space-time 
instantons describe the tunneling transitions between the states with differ- 
ent topological "winding" numbers v of the 577(2) <-> S 3 mapping [4, 5, 6, 7]. 
These topological numbers can be represented as the four-dimensional space- 
time integrals over the divergence of Chern-Simons current [8]. 

Topological fluctuations are believed to play an important role in the 
structure of QCD vacuum and in the properties of hadrons (for a review, see 
[9]). They also put in doubt the V and CV invariances of QCD ("the strong 
CV problem"). However until now all of the evidence for the topological 
effects in QCD from experiment, however convincing, has been indirect. In 
this paper I will summarize the arguments for the possibility to observe the 
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topological effects in QCD directly in the presence of very intense external 
electromagnetic fields. In particular, the coupling of topological gluon field 
configurations to electromagnetism induced by the axial anomaly leads to the 
separation of electric charges in the presence of a strong external magnetic 
field [10] ("the chiral magnetic effect" [11, 12, 13, 14]). The evidence for the 
chiral magnetic effect has been found recently from the numerical lattice QCD 
calculations [15]. The magnetic fields of the required strength can be created 
in heavy ion collisions [12, 16]. The direction of the produced magnetic 
field is orthogonal to the reaction plane of the collision. Very recently, STAR 
Collaboration presented an observation of the electric charge asymmetry with 
respect to reaction plane in relativistic heavy ion collisions at RHIC [17, 18]. 

2. The strong CP problem 

Let us begin with a brief introduction to the strong CP problem. Strong 
interactions within the Standard Model are described by Quantum Chromo- 
Dynamics, with the Lagrangian dictated by the SU (3) color gauge invariance: 



where G£ u and A afJi are the color field strength tensor and vector potential, 
respectively, g is the strong coupling constant, ipf are the quark fields of 
different flavors / with masses rrif, and t a the generators of the color £77(3) 
group in the fundamental representation. The Lagrangian (1) is symmetrical 
with respect to space parity V and charge conjugation parity C transforma- 
tions. 

However, these classical symmetries of QCD become questionable due to 
the interplay of quantum axial anomaly [19] and classical topologically non- 
trivial solutions - the instantons [3]. The axial anomaly arises due to the fact 
that the renormalization of the theory (1) cannot be performed in a chirally 
invariant way. As a result the flavor-singlet axial current = ipfy^ipf is 
no longer conserved even in the m — > limit: 



C = - 



\G^G a ^ v + V7 [*7 M (^ - igA afl t a ) - m f ] ip f , 



(1) 



<9 m J m5 = 2mfiil)fy 5 i/)f 




(2) 



where G atiV = \e^ vpa G apa ' . The last term in (2) is seemingly irrelevant since 
it can be written down as a full divergence, 



GfGV = d,K», (3) 
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of the (gauge-dependent) Chern-Simons current 



g 

Gapa — T^faPfAppA^v . (4) 



However this conclusion is premature due to the existence of instantons 
which induce a change in the value of the chiral charge Q 5 — J d 3 xK° asso- 
ciated with the topological current between t = — oo and t = +00: 

" = £~dt d ^ = 2N fq [G], (5) 

where 

q[G] = 32^ / d * xG «^ (6) 

is the topological charge; for a one-instanton solution, q — +1. 

In the presence of degenerate topological vacuum sectors, an expectation 
value of an observable O has to be evaluated by first computing an average 

(O) = [ DMD$]D[A]av(iSQCD)0(il,,j,A) (7) 

Jq 

over a sector with a fixed topological charge q, and then by summing over 
all sectors with the weight f(q) [20]. The additivity constraint 

/(5i + 52) = /(9i)/(g 2 ) (8) 
restricts the weight to the form 

f(q)=exp(i0q), (9) 

where 9 is a free parameter. Recalling an explicit expression (6) for q[G] one 
can see that this procedure is equivalent to adding to the QCD Lagrangian 
(1) Sqcd = J d A xC QC D a new term 

Ce = -^9 2 G^G afM/ . (10) 

Unless 9 is identically equal to zero, V and CP invariances of QCD are lost! 

Parity violation in strong interactions has been never detected, and strin- 
gent limits on the value of CP violating phase 9 < 3 x 10~ 10 follow from 
the experimental bounds on the electric dipole moment of the neutron [21]. 
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Perhaps the most appealing resolution of this puzzle is based on promoting 
the 9 parameter to a dynamical axion field [22, 23] emerging as a Nambu- 
Goldstone boson of an additional chiral symmetry [24]. We will see that 
even though in the physical vacuum 6 is equal to zero, the fluctuations of 
topological charge have an observable effect on QCD phenomena, and the 
concept of axions allows to formulate an appropriate effective theory. 

3. Chern-Simons diffusion 

At finite temperature, the transitions between the vacuum states with dif- 
ferent topological numbers can occur not only through quantum tunneling, 
but can also be induced by a classical thermal activation process, through 
a "sphaleron" [25]. In electroweak theory sphaleron transitions cause the 
baryon number violation and may be responsible for at least a part of the 
observed baryon asymmetry in the Universe [26]; for a review, see [27]. 
Sphalerons are also expected to play a role in QCD plasma [28] where they 
induce the quark chirality non- conservation. Unlike for the instanton transi- 
tions, the rate of the sphaleron transitions Y is not exponentially suppressed 
at weak coupling g, and in Yang-Mills theory with N colors is proportional 
to [29, 30, 31] 

r = const x (g 2 N)Hn(l/g 2 N) T 4 (11) 

(with a numerically large pre-factor [32]). Sphalerons describe a random 
walk in the topological number; in a volume V and after a (sufficiently long) 
time period T we get the topological number (z/ 2 ) = T V T. The diffusion of 
topological charge of course is expected to occur not only at weak coupling; 
while we cannot compute the corresponding rate analytically in QCD, lattice 
calculations indicate a large rate at experimentally accessible temperatures 
of T = 200^-300 MeV [32]. 

A valuable insight is offered also by the M = 4 maximally super-symmetric 
Yang-Mills theory (N=4 SYM) where the topological charge diffusion rate 
can be evaluated in the strong coupling limit through the AdS/CFT corre- 
spondence [33, 34, 35], with the following result [36]: 

T = (# 2 A07(2567r 3 ) T 4 , (12) 

which shows that the topological transitions become more frequent at strong 
coupling, even though the dependence on the coupling is weaker than sug- 
gested by (11). Note that the large N behavior is the same in the weak and 
strong coupling limits (~ N°). 
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An explicit example of a topological solution in weakly coupled AdS 5 x S 5 
supergravity (that is dual to the large N, strongly coupled N=4 SYM) is 
provided by the D-instanton [37, 38] that can be obtained from solitonic D- 
brane solutions by wrapping them around an appropriate compact manifold. 
D-instanton can be viewed in the string frame as an Einstein-Rosen worm- 
hole connecting two asymptotically Euclidean regions of space-time, with 
the Ramond-Ramond (R®R) charge flowing down the throat of the worm- 
hole [37]. It describes a violation of the conservation of a global charge in 
physical processes [37]. Since in the AdS/CFT dictionary the R®R scalar 
(the axion) of supergravity is dual to the 8 angle of N=4 SYM field theory 
[39], this flow of R®R charge down the throat of the wormhole describes the 
change of topological charge in the field theory description. D-instantons 
have recently been considered [40] as a source of multiparticle production in 
high energy collisions in strongly coupled N=4 SYM. 



4. Topologically induced effects in electrodynamics: 
Maxwell-Chern-Simons theory 

4-1. The Lagrangian 

Let us begin by coupling the theory (1) to electromagnetism; the resulting 
theory possesses 577(3) x U{\) gauge symmetry: 

£qcd+qed = -^G^G aflu + ^2i>f [«7 m («9m - igA a ^t a - iq f A^) - m f ] tp f - 

f 

-^g^G^-^F^, (13) 

where and are the electromagnetic vector potential and the corre- 
sponding field strength tensor, and qf are the electric charges of the quarks. 

Let us discuss the electromagnetic sector of the theory (13). Electromag- 
netic fields will couple to the electromagnetic currents J M = X)/ 9/^/7/* V*/- 
In addition, the term (10) will induce through the quark loop the coupling of 
FF to the QCD topological charge (see Fig.l). We will introduce an effective 
pseudo-scalar field 9 = 9(x,t) (playing the role of the axion field) and write 
down the resulting effective Lagrangian as 

Cues = -\f<"F^ - A^.P - | 9F^F^, (14) 
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where 

c = ^gje 2 /(27r 2 ). (15) 
/ 

This is the Lagrangian of Maxwell-Chern-Simons, or axion, electrody- 
namics that has been introduced previously in [41, 42, 43]. If 9 is a constant, 
then the last term in (14) represents a full divergence 



F^F, 



pV 



°P J CS 



(16) 



of the Chern-Simons current 



TP _ A rp 

°CS — f^-v 1 pat 



(17) 



which is the Abelian counterpart of (4). Being a full divergence, this term 
does not affect the equations of motion and thus does not affect the electro- 
dynamics of charges. 

The situation is different if the field 9 = 9(x,t) 
varies in space-time. Indeed, in this case we have 



9F^F, 



flU 



0d P J CS 



dp [er cs \ - dpOJ£s- (is) 



The first term on r.h.s. is again a full derivative and 
can be omitted; introducing notation 



P^ = d fl 9 = (M,P) 



(19) 



we can re- write the Lagrangian (14) in the following 
form: 



1 



MCS 



ppyp 



pv 



P 7 M 



(20) 




Since 9 is a pseudo-scalar field, P^ is a pseudo- 
vector; as is clear from (20), it plays a role of the po- 
tential coupling to the Chern-Simons current (17). 
However, unlike the vector potential A^, P^ is not 
a dynamical variable and is a pseudo-vector that is 
fixed by the dynamics of chiral charge 



m our case, 



Figure 1: The quarks cou- 
ple the topological Chern- 
Simons current to the elec- 
tromagnetic field through 
the axial anomaly. 



determined by the fluctuations of topological charge in QCD. 

In (3 + 1) space-time dimensions, the pseudo- vector P M selects a direction 
in space-time and thus breaks the Lorentz and rotational invariance [42] : the 
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temporal component M breaks the invariance w.r.t. Lorentz boosts, while 

— * 

the spatial component P picks a certain direction in space. On the other 
hand, in (2 + 1) dimensions there is no need for the spatial component P 
since the Chern-Simons current (17) in this case reduces to the pseudo-scalar 
quantity e vpa A u F pa , so the last term in (20) takes the form 

AC = c Mt UfX7 A u F pa . (21) 

This term is Lorentz-invariant although it still breaks parity. In other words, 
in (2 + 1) dimensions the vector P can be chosen as a 3- vector pointing in the 
direction of an "extra dimension" orthogonal to the plane of the two spatial 
dimensions. When added to the Maxwell action, (21) generates a mass of 
the photon which thus becomes "topologically massive". This illustrates an 
important difference between the roles played by Chern-Simons term in even 
and odd number of space-time dimensions. 

4-2. Maxwell- Chern-Simons equations 

Let us write down the Euler-Lagrange equations of motion that follow 
from the Lagrangian (20), (17) (Maxwell-Chern-Simons equations): 

d p F^ = J v - P P F^. (22) 

The first pair of Maxwell equations (which is a consequence of the fact that 
the fields are expressed through the vector potential) is not modified: 

dpF* 1 " = r. (23) 

— * 

It is convenient to write down these equations also in terms of the electric E 
and magnetic B fields: 

V x B - ^ = J + c (MB -Pxl), (24) 
V • E = p + cP ■ B, (25) 

VxE + — =0, (26) 

V • B = 0, (27) 

where (p, J) are the electric charge and current densities. One can see that 
the presence of Chern-Simons term leads to essential modifications of the 
Maxwell theory. We will now examine a few examples (some old and some 
new) illustrating the dynamics contained in Eqs(24)-(27). 
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4.3. The Witten effect 




Figure 2: Magnetic monopole at finite 9 angle acquires an electric charge ~ e8/ir that is 
localized on the domain wall where the value of 9 changes from zero in the core of the 
monopole to some value 9 ^ away from the monopole (the domain wall is shown by the 
gray ring) - the Witten effect [44]. 

Our first example will be the celebrated Witten effect [44]: magnetic 
monopoles at finite 9 angle acquire electric charge and become "dyons". 
Our treatment of this effect will follow Wilczek [41]. Consider a magnetic 
monopole in the presence of a finite 9 angle. In the core of the monopole 
9 = 0, and away from the monopole 9 acquires a finite non-zero value - there- 
fore within a finite domain wall we have a non-zero P = V# pointing radially 
outwards from the monopole (see Fig. 2). According to (25), the domain wall 
thus acquires a non-zero charge density cV# ■ B. An integral along P (across 
the domain wall) yields J dl 89/ dl = 9, and the integral over all directions of 

— * 

P yields the total magnetic flux $. By Gauss theorem, the flux is equal to 
the magnetic charge of the monopole g, and the total electric charge of the 
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configuration is equal to 



, = c,, = s'' = S5 ,w = e S' (28) 

where we have used an explicit expression (15) for the coupling constant c, 
as well as the Dirac condition ge = An x integer. 



5. Chiral magnetic effect 

5.1. Charge separation 

Consider now a configuration shown on Fig. 3 where an external magnetic 
field B pierces a domain with 9 ^ inside; outside 9 = 0. Let us assume first 
that the field 9 is static, 9 = 0. Assuming that the field B is perpendicular 
to the domain wall, we find from (25) that the upper domain wall acquires 
the charge density per unit S of [11] 

= + c 9B (29) 

up 

while the lower domain wall acquires the same in magnitude but opposite in 
sign charge density 

(|) = - c 9B (30) 

\ / down 

Assuming that the domain walls are thin compared to the distance L between 
them, we find that the system possesses an electric dipole moment 

d e = c9(B.S)L = ^aj (e£) (^) L; (31) 

in what follows we will for the brevity of notations put qj — 1; it is easy 
to restore this factor in front of e 2 when needed. 

Static electric dipole moment is a signature of V, T and CP violation (we 
assume that CVT invariance holds). The spatial separation of charge will 
induce the corresponding electric field E = c 9 B. The mixing of pseudo- 
vector magnetic field B and the vector electric field E signals violation of V, 
T and CP invariances. 

The formula (31) allows a simple interpretation: since eB /2ix is the trans- 
verse density of Landau levels of charged fermions in magnetic field B, the 
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Figure 3: Charge separation effect - domain walls that separate the region of 9 ^ from 
the outside vacuum with 6 = become charged in the presence of an external magnetic 
field, with the surface charge density ~ ed/ir ■ eB/2ir. This induces an electric dipolc 
moment signaling V and CV violation. 

floor of the quantity eB ■ S/2tt (i.e. the largest integer that is smaller than 
eB ■ S/2n) is an integer number of fermions localized on the domain wall. 
Each fermion species contributes independently to this number as reflected 
by the factor Nf. Again we see that the electric dipole moment (31) arises 
from the electric charge q ~ e8/7r that is induced on the domain walls due 
to the gradient of the pseudo-scalar field 6. 

If the domain is due to the fluctuation of topological charge in QCD 
vacuum, its size is on the order of QCD scale, L ~ Aqq D , S ~ ^qcd- This 
means that to observe an electric dipole moment in experiment we need an 
extremely strong magnetic field eB ~ Aq CD . Fortunately, such fields exist 
during the early moments of a relativistic heavy ion collision [12]. Here 
we have assumed that the domain is static; this approximation requires the 
characteristic time of topological charge fluctuation r ~ 1/9 be large on the 
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time scale at which the magnetic field B varies. This assumption is only 
marginally satisfied in heavy ion collisions, and so we now need to consider 
also the case of 9 ^ 0. 

5.2. The chiral induction 

— * 

Consider now the domain where \P\ <C 9, i.e. the spatial dependence of 
9(t, x) is much slower than the dependence on time [12], see Fig. 4. Again, we 
will expose the domain to an external magnetic field B with V x B = 0, and 
assume that no external electric field is present. In this case we immediately 
get from (24) that there is an induced current [13] 

J = -cMB = --^-9B. (32) 

— * 

Note that this current directed along the magnetic field B represents a V— , 
T— and CP— odd phenomenon: the electric current is a vector while the 
magnetic field is a pseudo- vector. This current is of course absent in the "or- 
dinary" Maxwell equations that respect V and CP invariances. Integrating 
the current density over time (assuming that the field B is static) we find 
that when 9 changes from zero to some 9 ^ 0, this results in a separation of 
charge and the electric dipole moment (31). 

Let us discuss the meaning of formula (32) in more detail. To do this, let 
us consider the work done by the electric current; to obtain the work per unit 
time - the power P - we multiply both sides of (32) by the (static) electric 

— * 

field E and integrate them over the volume (as before, we assume that 9 does 
not depend on spatial coordinates): 

p = J d 3 x J ■ E = -9 ^ J d 3 x E ■ B = -9 Q 5 , (33) 



where 



e 2 



Q 5 = — 2 I dt d 6 x E ■ B (34) 



is the chiral charge. The meaning of the quantity on the r.h.s. of (33) can 
be revealed with the help of the following well-known quantum-mechanical 
analogy. The vacuum wave function 

|0) = J]exp(i0Q5) IQs) (35) 

Qr, 
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B ♦ 




Figure 4: The chiral magnetic effect - inside a domain with 9 ^ an external magnetic 
field induces an electric current J ~ eO/ir ■ eB/2ir. 9 =^ indicates the change of the chiral 
charge inducing an asymmetry between the left- and right- handed fermions inside the 
domain. Note that the current J ~ B is absent in Maxwell electrodynamics. 

is analogous to the Bloch wave function of electron in a crystal, with 9 playing 
the role of electron's quasi-momentum, and - the role of coordinates of 
atoms in the crystal. The derivative of the "momentum" 9 thus plays the role 
of the force that we assume is constant in time and Q5 - of the dimensionless 
distance; Q 5 is thus the velocity. The formula (33) is therefore simply the 
classical expression 

Power = Force x Velocity, 
with the force acting along the "extra dimension" of the chiral charge Q$. 

Let us stress that we have derived (32) from the Maxwell-Chern-Simons 
equations in the absence of external charges and currents. Therefore, the 
current (32) should arise from the decay of the vacuum state induced by 
the change of chirality in the presence of magnetic field. Since the vacuum is 
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Left 

V ) 



Right 



Figure 5: Dirac cones of the left and right fermions. In the presence of the changing chiral 
charge there is an asymmetry between the Fermi energies of left and right fermions hl 
and fi R : [i L - (1r = 2^ 5 = 29. 

supposed to have zero energy, what powers the "chiral magnetic" current? To 
answer this question, consider the physical vacuum of chiral Dirac fermions: 
all negative energy states of both left and right fermions are occupied, and 
all positive energy states are empty. The anomaly relation 

Q 5 = N R - N L (36) 

tells us that in the presence of chiral charge there should be an asymmetry 
between the numbers of left and right fermions, and thus an asymmetry in 
the corresponding Fermi momenta pf,R and Fermi energies Hl.r = c Pl r, see 
Fig. 5.2. 

Let us compute the densities of left and right fermions in (anti-) parallel 
electric and magnetic fields; for definiteness let us choose the negative sign 
for the product E ■ B < 0. The presence of magnetic field B aligns the spins 
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of the positive (negative) fermions in the direction parallel (anti-parallel) to 

— * 

B. In the electric field E the positive fermions will experience the force eE 
and will move along E; therefore their spin will have a positive projection on 
momentum, and we are dealing with the right fermions. Likewise, the nega- 
tive fermions will be left-handed. After time t, the positive (right) fermions 
will increase their Fermi momentum to p 1 ^ = eEt, and the negative (left) will 
have their Fermi momentum decreased to p£ = — p^. The one-dimensional 
density of states along the axis z that we choose parallel to the direction of 
fields E and B is given by dN R /dz = p R /2n. In the transverse direction, the 
motion of fermions is quantized as they populate Landau levels in the mag- 
netic field. The transverse density of Landau levels is d 2 N R /dxdy = eB/2n. 
Therefore the density of right fermions increases per unit time as 

d ' NR - e " f n cvr\ 
dTdv-J^y E - B - (37) 

The density of left fermions decreases with the same rate, d 4 Ni/dt dV = 
—d 4 N R /dt dV. The rate of chirality Q 5 = N R — N L generation is thus 

d A Q^ e 2 - - . . 

ww = ^ E - B - < 38 > 

which is exactly the 4D density of topological charge as given by (34). 

Transferring a fermion from left Fermi surface to the right one costs an 
amount of energy equal to fi R — fii, see Fig. 5.2. The power P (work per 
unit time) of the current [45, 13] is thus equal to 

p = J d 3 x J-E=(hr-hl) j d 3 x E-B (39) 

This coincides with (33) if we identify 

10 — hl fi R = 2fi 5 , (40) 

where we have defined a chiral chemical potential /is. Again, this matches 
our interpretation of 6 as of the "force" that creates the net momentum 
115/c along the "extra dimension" of topological charge Q5. This anomalous 
mechanism of current generation is analogous to the one invoked by Witten 
[46] to describe the particle acceleration in cosmic strings. Related condensed 
matter examples have been considered in [47, 48]. 
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5.3. Charge separation at finite baryon density and vorticity: 
the chiral vortical effect 
At finite baryon density, the charge separation and the chiral induction 
can occur even without an external magnetic field if the angular momentum 
is present [11, 10]. Indeed, let us introduce a matter velocity field V^; then 
at finite baryon density \i Eq. (20) acquires the following additional term: 

Cb = -N c ^ ■ e^d^e {d x V v )A a] (41) 

this term has been used in the studies of the axial current in cold dense 
quark matter [49, 50]. It is easy to see that it also induces the electric charge 
density on the topological domain walls [11]: 

PB = ^ = ■ (d 3 V k ) = N c ^ 2 ■ (W ■ n) , (42) 

as well as an electric current directed along the along the angular velocity Q of 
the rotating system, 2eij k VL k = (diVj — djVi). In a recent paper [51] it has been 
shown how this phenomenon manifests itself in relativistic hydrodynamics. 
Since the plasma produced in off-central heavy ion collisions should possess 
vorticity, this "chiral vortical effect" may be an important source of charge 
separation, especially at relatively low energies when the net baryon density 
of the produced matter is significant. 



6. Experimental observations 



As we discussed above, the fluctuations of topological charge in the pres- 
ence of strong magnetic field should induce the separation of electric charges. 
In heavy ion collisions, the magnetic fields of required strength are produced 
naturally by the electrically charged ions in the initial state, spectators in 
the final state, and due to the electric charge asymmetry in the distributions 
of the produced hadrons [12, 16]. The produced magnetic fields are oriented 
perpendicular to the reaction plane (along the system's angular momentum); 
therefore the chiral magnetic effect under discussion should result in the sep- 
aration of the electric charge with respect to the reaction plane [10]. 

The azimuthal distribution of charged hadrons produced in heavy ion 
collisions can be expanded in Fourier harmonics in the following way: 

~ 1 + 2vi cos(A0) + 2v 2 cos(2A0) + ... + 2a± sin(A0) + (43) 
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Figure 6: The STAR Collaboration result on the charge-dependent azimuthal correlations 
in Au-Au collisions at y/s = 200 GeV per nucleon pair at RHIC; from [17, 18]. Horizontal 
axis is the centrality of the collision (the fraction of the inclusive inelastic cross section); 
the average impact parameter decreases towards the right, as does the magnetic field. Also 
shown are the predictions of various Monte-Carlo models of heavy ion collisions. 



where A0 = — ^rp is the angle with respect to the reaction plane - the 
plane which contains the impact parameter and beam momenta. Note that a 
typical relativistic heavy ion collision produces several thousand hadrons, so 
the reaction plane can be reliably identified in each event. The coefficients v± 
and v 2 measure the strength of so-called directed and elliptic flow. The chiral 
magnetic effect causes the hadrons with opposite electric charges to be pro- 
duced preferentially on opposite sides of the reaction plane; this corresponds 
to non-zero coefficients a + and a_ that should have the opposite signs; this 
form of the angular distribution follows from the structure of P-violating 
amplitude [10]. 

The number of charged hadron tracks in the single event (although suffi- 
cient to determine the reaction plane) is not large enough to allow a statis- 
tically sound extraction of the coefficients a±, so one has to sum over many 
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events. However since there is no global violation of V and CP invariances in 
QCD, the sign of the charge asymmetry should fluctuate from event to event 
and so when averaged over many events, (a+) = (a+) = 0. The way out of 
this dilemma has been proposed by Voloshin [52] who suggested to extract the 
cumulant {a a ap) by measuring the expectation value of (sin(A0 a ) sin(A<^g)). 
The variable that was proposed in [52]: 

(cos(0 Q + 4>f3 — 2^ RP )} = (cos A<fi a cos Atfip) — (sin A0 Q sin A<^g) (44) 

has an added benefit of not being sensitive to the reaction plane-independent 
backgrounds that cancel out in (44). The variable (44) can be measured with 
a very high precision, and is directly sensitive to the parity-odd fluctuations; 
the price to pay however is that the observable itself is parity-even, and so 
one has to carefully examine all possible physics backgrounds. 

The preliminary results have been reported by STAR Collaboration in 
[53, 54]. Very recently, the STAR Collaboration has presented a conclusive 
measurement of (44) that amounts to the observation of charge-dependent 
azimuthal asymmetries in heavy ion collisions at RHIC [17, 18]; one of the re- 
sults is shown in Fig.6. One can see that the same-charge and opposite-charge 
cumulants (44) differ in a way that is very significant statistically. Numerous 
mundane backgrounds have been examined, and none of them could explain 
the observed effect so far [17, 18]. The predictions of various Monte Carlo 
models of heavy ion collisions are also shown in Fig.6; these models (while 
successful in reproducing the global features of heavy ion collisions) fail in 
explaining the observed effect. Needless to say, one has to continue looking 
for alternative explanations of the STAR result, and to extend the experi- 
mental study of charge asymmetries; however in any case the observed effect 
is clearly novel and very intriguing. 

7. Summary and discussion 

In this paper we discussed how the fluctuations of topological charge in 
the presence of strong magnetic field induce the separation of electric charge 
along the axis of magnetic field (the chiral magnetic effect) [10, 11, 12, 13, 
14]. The evidence for the chiral magnetic effect from lattice QCD has been 
found recently in [15]. The behavior of the effect at strong coupling in the 
Sakai-Sugimoto model and related theories has been explored through the 
AdS/CFT correspondence in [55, 56, 57]. Recently, the STAR Collabora- 
tion at RHIC has presented an observation of charge-dependent azimuthal 
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asymmetries [17, 18]. None of the mundane backgrounds could explain the 
observed effect so far [17, 18]. It is clear that a dedicated experimental pro- 
gram is necessary to understand fully this intriguing observation. In particu- 
lar, one could extend the present studies of charge asymmetries to identified 
hadrons and to look for the possible weakening of the observed effect at low 
energies where the energy density is under the threshold for forming the QCD 
plasma. Both studies have already been planned at RHIC. 

Here we discussed the sphaleron transitions as a way of generating topo- 
logical charge in the QCD plasma. There exists however a number of other 
possibilities: the metastable parity-odd domains [58, 59, 60, 61] that emerge 
in the large N description of plasma in QCD and related theories; and the 
dyons that have been identified in the plasma in lattice QCD [62, 63], to 
mention a few. Whatever the specific "microscopic" origin of the topological 
fluctuations, the fact of their existence in the plasma is firmly established by 
now in numerical lattice simulations. The recent study [64] finds that the 
topological charge density in QCD plasma at temperatures T c < T < 2T C is 
correlated over distances of order 1/T, and that these correlations are much 
stronger than in the vacuum; for a review of the earlier work, see [65] . It is im- 
portant to note that topological fluctuations exist already at the early stages 
of the collision prior to thermalization, both in the weak [66, 67, 68, 69, 70] 
and strong [40] coupling descriptions. 

The mechanism considered here requires a sufficiently large energy density 
for the topological transitions to turn on, and for the quarks to separate 
by a distance comparable to the system size - therefore, there has to be a 
deconfined phase. In addition the system has to be in a chirally symmetric 
phase - in a chirally broken phase, the chiralities of quarks could flip easily 
causing dissipation of the induced current. The experimental and theoretical 
studies of parity-odd charge asymmetries in heavy ion collisions can thus 
bring us closer to the detection of a change in the physical properties of the 
vacuum induced by creating a high energy density over an extended volume 
[71]. 

Note added: Very recently, Blum and collaborators [72] have reported 
on the first study of the chiral magnetic effect in lattice gauge theory with 
dynamical 2 + 1 flavors of quarks (in the domain wall formulation). 
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